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I. HIST ORY AND INTRODUCTION 



The U. S. Navy is a hardware oriented organization 
reguiring a wide variety of spare and replacement parts. 
The Naval Supply Systems Command (NAVSUP) is the organiza- 
tion within the Navy responsible for maintaining a suffi- 
cient stock of parts and replenishing this stock when 
quantities run low. NAVSDP is keenly interested in the 
guality of the models it uses to provide support for the 
Navy’s weapons systems. Faculty of the Naval Postgraduate 
School were asked to examine the existing wholesale and 
retail provisioning and replenishment models and to provide 
comments on their correctness and suggestions for possible 
improvements. In response, a report by Richards and 
McMasters [Ref. 1] addresses the problems and various solu- 
tion techniques. This thesis examines a specific solution 
technique that may be applicable to one of the problems. 
The following is a brief synopsis of the existing models and 
the environment in which they operate. 

Within NAVSUP there are two centers that use provi- 
sioning models, 1) Ships Parts Control Center (SPCC) and 2) 
Aviation Supply Office (ASO) . Presently the Inventory 
Control Points (ICPs) within the centers have two retail 
provisioning models which generate, 1) COSALS (Coordinated 
Shipboard Allowance Lists) and 2) AVCALS (Aviation 
Coordinated Allowance Lists). These models supply suffi- 
cient parts support for a short term ship deployment and try 
to provide protection from replacement part depletion until 
the wholesale replerishmen t models can buy or repair the 
needed spare parts. 

Weapon system procurement requires that provisioning 
support be available by a set preliminary operational 
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capability (POC) date. The initial support requirements are 
determined by a wholesale provisioning model. The support 
must be comprehensive enough to supply a new weapon system 
until the replenishment buy can be made and delivered. It 
should be noted that the provisioning lead time is usually 
two to two and one half years. Variability in production 
and installation causes the quality of the forecasted demand 
to be poor. Hence, the provisioning buy may result in large 
quantities of excess parts in the supply system, or there 
may be many stockouts before the replenishment buy becomes 
available. The Assistant Secretary of Defense 

(Installations and Logistics) has provided guidance [Ref. 2] 
to help prevent large excesses. From this instruction a 
third model (COSDIF) was developed to conservatively esti- 
mate the quantities of new parts required for wholesale 
provisioning. For existing equipment no additional stock is 
procured for the new weapons system. If the increase in 
demand from the new system is large enough, the inventory 
manager can anticipate the increase and start a replenish- 
ment buy in sufficient time to meet the need. Otherwise, 
inventory optimization models may be used to control the 
stock. Wholesale replenishment for a new system is initi- 

ated at different times by the two centers within NAVSUP. 
For either center the reorder point for replacement parts is 
subject to change as a function of the actual demand the 
weapon system places on the supply system. The reorder 
point for an existing system occurs at part depletion levels 
set by NAVSOP. 

Budgeting constraints for support of existing systems 
are partially driven by the peculiarities of the Department 
of Defense Planning, Programming and Budgeting System 
(PPBS). Those purchases for new weapon systems are funded 
as part of the procurement process. For existing systems, 
each Service requests from Congress sufficient funding to 
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meet specific levels of readiness and availability set by 
the respective Service Head. NAVSUP receives this direction 
from the Chief of Naval Operations (CNO) and then priori- 
tizes the parts required with a procedure called "Variable 
Threshold". This technique considers several inputs 
peculiar to the weapon system to develop the priority list 
for items to be purchased. The quantity (depth) to be 
purchased is determimed from historical data or from engi- 
neering estimates of the number of anticipated failures for 
a period of time equal to one procurement lead time plus cne 
quarter. Proceeding from the top of the priority list the 
parts are purchased to their ’depth’ until the budget is 
expended. Those parts unfunded for purchase are placed in a 
file for Unfunded Requirements. 

Obviously, the optimal quantities to purchase, as a 
function of the inputs, and the maximum utilization cf the 
budget are of high interest. Howard in £Ref. 3] addresses 
several different approaches that provide solutions for this 
type problem. The thrust of this thesis is a variation of 
one of those solution techniques for a similarly defined 
provisioning problem. The objective is to minimize the Mean 
Supply Response Time (MSRT) subject to a budget constraint. 

The overall MSRT is a nonlinear composite function of 
the purchase quantity for all items. The individual item's 
MSRT functions are also nonlinear with performance as a 
function of four variables peculiar to the item type. The 
budget constraint is a linear function with constant costs 
for each item and a set budget limit. There are no cost 
advantages for large magnitude buys or item combination 
buys. The individual costs and return functions of the MSRT 
lend themselves to formulation as a dynamic programming 
(DP) problem. However, the problem is large and complex for 
a standard recursive dynamic approach due to the required 
number of stages, the magnitude of the budget and the range 
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a variation of the normal D? 



of costs. Consequently, 
approach was developed that significantly reduces the 
required computations. This variation shall be referred to 
as 'the funnel' because of its funneling characteristics for 
computation of the resource variable bounds. Further expla- 
nation of its construction and procedure may be found in 
Chap 2&3. Modifications of an existing DP computer program 
described in [Bef. 3] to operate as the DP variation makes 
the 'funnel' method a usable approach considering speed and 
ease of manipulation. Additionally, the DP variation has 
the advantage of providing an optimal integer solution. In 
this report 'the funnel' will be compared to a full DP 
procedure and a computerized marginal analysis approach 
which operates very efficiently but which may not generate 
the optimal solution. 
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II. SOLUTION TECaNIQOES 



This chapter will discuss three methods for solving the 
provisioning problem described in Chapter 1, 1) marginal 

analysis, 2) pure dynamic programming, and 3) a modified 
dynamic programming procedure. The first technique provides 
a quick solution that is a good estimate but may not be 
optimal. The second technique provides a global optimal 
solution but requires a large amount of computer time for 
any moderate sized problem. The third technique provides a 
local, and possibly a global, optimal solution and requires 
significantly less time than option two. 

The following definitions will be used in the formula- 
tion and discussions in this chapter. Capital letters 
represent the variable while lower case letters indicate a 
subscript. In the expression Xi-1 the quantity i-1 is the 
subscript although it appears on the same line as the symbol 
X. 



N 

3 

ci 

Ei 

Ki 



Li 

Si 



Ti 

P'i (Si) = 



the total number of items considered for 
provisioning 

pre-assigned maximum budget amount 
the unit cost of item i 
the essentiality code for item type i 
the fixed time required to satisfy a 
requisition for item i if stock is 
available 

the demand rate for item i (lambda) 
the decision variable for the number 
of items of type i 

the procurement lead time for item i 
the probability of Si demands for item i 
during the provisioning interval 
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Pi (Si) = the probability of Si or fewer demands 
for item i daring the provisioning 
in terval 

Zi (Si) = the performance measure for item i when 
Si units are stocked. 

The objective of this problem is to minimize the Mean 
Supply Response Time (MSRT) as a function of Si. MSRT is a 
widely used measure cf supply performance because of its 
role in determining availability and because it is an indi- 
cator of the success of a supply system in meeting response 
time goals. It considers two major factors in computation, 

1) the likelihood of satisfying demands from stock on 
hand and 

2) the length of the delay in satisfying demands 
when the system runs out of stock. 

It may be represented by the following expression, 

MSRT=min £Ei*Li*Ti*Zi (Si) /£Ei*Li* Ti (2.1) 

i = 1, ,N, 

where the likelihood adjustments for conditions one and two 
are factors in the Zi(Si) function. This function repre- 
sents the performance measure of item i when Si items are 
considered and is defined by the following equation for the 
case where demands during the provisioning leadtime are 
assumed to be Poisson distributed: 

Zi (Si) =Ki + ( (1-Pi (Si) ) * (Li *Ti-2 Si+ Si * { Si+ 1) (2.2) 

/Li*Ti) +P’ i (Si) (Li*Ti- Si) ) /2Li) . 

The constraint for this problem is the budget or the 
resource in dollars. It is linear and may be represented by 
the following inequality. 
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£ ci*Si<=B 



i= 1 



i • • • f 



(2.3) 



N. 



The following subsections present solution technigues which 
are viable procedures for solving the above formulated 
provisioning problem. 

A. MARGINAL ANALYSIS 

The marginal analysis method is a guick, simple and 
effective approach for generating a solution to the previ- 
ously defined provisioning problem. For the problem consid- 
ered here, the marginal analysis approach generates 
solutions which may not be optimal. However, the solutions 
generated by the procedure do possess some very desirable 
properties which are stated below. 

Let Zi (Si) be the essentiality weighted MSRT for item i 
when Si units are stocked. Let ^Zi (Si) =Zi (Si-1) -Zi (Si) be 
the reduction in essentiality weighted MSRT achieved by 
increasing the stockage level for item i from Si-1 units to 
Si units. Let S 1 = ( S1 # S2, . . . , SN) be the vector of alloca- 
tions for the N items; i.e. Si is the number of units allo- 
cated to item i. The marginal approach begins with a budget 
of B dollars and an allocation vector S’ = (0 , 0, . . . , 0) . It 
computes £Zi(1)/Cl for i=1,2,...,N and selects that item for 
which this marginal benefit to cost ratio is greatest. 
Without loss of generality , assume that the item selected 

is item 1. This requires an expenditure of Cl dollars and 

results in the allocation vector S 1 = ( 1 , 0, 0 , 0) . The 

procedure next computes £Z 1 (2) /Cl, the marginal benefit to 
cost ratio achieved by increasing the stockage for item 1 by 
one unit from a level of one to two units. The maximum 
value is then selected from among the ratios in 2.4. 

An additional unit of the item, say item j, for which the 
maximum occurs is selected and the remaining budget is 
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(2.4) 



Az 1 ( 2 ) , Az 2 ( 1 ) , Az 3(1) .... , AZN( 1 ) 

Cl C 2 C 3 CN 

decremented by C j dollars. The process continues in this 
manner allocating a single unit iteratively until the budget 
is expended or no additional units can be purchased. At 
each step, the ratios compared are 

Az 1 (S 1 + 1) , Az 2 (S 2 + 1 ) , ... , AzN ( SN+1) (2.5) 

Cl C 2 CN , 

where the current allocation vector is S’ = (S 1 , S2 , . . .SN) . 
Since only one of the ratios changes at each step of the 
process, the computations required are very few. 

It can be shown [Eef. 4], that for the problem consid- 
ered here the solution is undominated. That is, if B(S’) is 
the budget required to fund the allocation vector S’, and if 
S’ 1 is any other allocation vector such that 3 (S * * ) <=B (S’ ) 
then the overall essentiality weighted MSRT for S’ is 
smaller than that for S ,f . If it should happen that B(S') = 
B, the original budget, for some step of the marginal anal- 
ysis procedure, then the solution S' is optimal. 

In actual practice, what usually happens in the marginal 
analysis procedure is that at some step purchase of the item 
with the largest marginal benefit to cost ratio requires an 
expenditure of more money than that which remains. The 
allocation vector at that point is not feasible and various 
heuristics are usually applied to continue the process. The 
most popular heuristic is to backtrack to the previous 
feasible solution and select that item which has the maximum 
benefit to cost ratio among those items with unit costs no 
larger than the budget remaining. On applying such heuris- 
tics, one can not claim that the final solution so obtained 
is undominated. However, one can make use of the undomi- 
nated property to obtain a bound on the optimal value of the 
objective function. 
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Let Z (S S) be the objective function value at the occa- 
sion at which the marginal analysis procedure first gener- 
ates a solution for which the budget required is equal to or 
greater than B. If equality exists, SS is optimal. Suppose 
B (S S) =BS>B and let S* (X) be the optimal solution for a 
budget of X dollars. Then, since SS is undominated, 
Z (SS) <=Z (S* (BS) ) <= Z (S* (B) ) , which is the value of the objec- 
tive function at the optimal solution. Thus Z(SS) is a 
lower bound on the optimal value of the objective function. 
One can thus assess the potential benefit of a more exten- 
sive search for the optimal solution by comparing the actual 
performance obtained at any stage with the lower bound 
discussed above. 

A computer program was written by Richards to implement 
this method. A sample of this Marginal Analysis program 
output in table I shows the iteration number, which items 
are purchased, the numbers of items purchased, the benefit 

to cost ratios and the budget remaining. The table shows 
the sequence in which the items are purchased. Note the 
decreasing benefit to cost ratio. The program selects item 
number one five times in a row before the first unit of item 
two provides a better marginal return. The last attempt to 
purchase chooses an item that drives the remaining budget 
negative. The program automatically re-evaluates the 
marginal returns and then chooses the next largest benefit 
to cost ratio which in this case drives the remaining budget 
to zero. As previously discussed this solution may not be 
opti mal. 

The Marginal Analysis program has several other capabil- 
ities based on the same approach. There is a job selection 
menu that provides a choice of 1) spares determination, 2) 
budget determination or 3) kit evaluation. The first option 
selects the quantities to be purchased within a given 
budget. The second option solves a related problem by 
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TABLE I 

Sample Output from the Marginal Analysis Program 






*********************************************** 

* METHOD: FIXED EUDGET = 205-00 * 

*********************************************** 

*********************************************** 



STEP ITEM SI BATIO BB 



1 


1 


1 


0.480808496 


200.00 


2 


1 


2 


0.365659833 


195-00 


3 


1 


3 


0-260618091 


190.00 


4 


1 


4 


0- 172421157 


185.00 


5 


1 


5 


0. 105280340 


180-00 


6 


2 


1 


0.0735758543 


170.00 


7 


3 


1 


0.0600002967 


155.00 


8 


1 


6 


0.0591956675 


150.00 


9 


3 


2 


0. 0533369593 


135.00 


JO 


3 


3 


0-046683754 1 


120.00 


11 


3 


4 


0- 04C0909968 


105. 00 


12 


3 


5 


0.0336193480 


90.00 


13 


1 


7 


0.0306577347 


85.00 


14 


3 


6 


0.0273999237 


70.00 


15 


3 


7 


0.0216008648 


55.00 


16 


2 


2 


0-0207276568 


45.00 


17 


3 


8 


0. 0164023377 


30.00 


18 


1 


8 


0-0146531016 


25.00 


19 


3 


9 


0.0119544677 


10.00 


20 


3 


10 


0.0083406679 


-5.00 


20 


1 


9 


0.0064818673 


5.00 


21 


1 


10 


0.0026624922 


0.0 



****************************************** 


MINIMIZE 


MSBT 


EXCESS: 0.0 


I 


SI 


MSRI 


1 


10 


0. 190 


2 


2 


10.396 


3 


9 


12. 176 



OVEBALL 5-185 

REALIZED PERFORMANCE = 5.3852 

BOUNDS { 6.0033, 6.3120) 
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B. STANDARD DYNAMIC PROGRAMMING 



The standard dynamic programming approach can be applied 
to this type of provisioning problem. The following sample 
problem illustrates how the concept may be applied. A 
single stage problem corresponding to a single item may be 
characterized as in fig 2.1. 




Figure 2.1 Example of a One-Stage Decision Problem. 



In figure 2.1, 



X is the input state variable 
available for use at this 



which is the 
stage. 



resource 



S is the decision variable that represents the decision 
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within the stage. 

E is the stage return function that produces a return 
fcased on the input X and the decision S. 

Y is the output state which is a function of X and S. 

T is the stage transformation function that expresses 
the components of the output state Y as a function 
of the input state and decisions, Y=T(X,S). 

An after an Xi or Si variable indicates the optimal 

value for that stage, [Eef. 5 p2 2-23.]. 

The one-stage initial-state optimization problem is to find 
the best stage return as a function of the input state X. 
The optimal return from the stage will be denoted as F (X) 
and the optimal decision as S* (X) . Thus, 



F (X) = min E (X,S (X) ) = R(X,S* (X) ) , (2.6) 

S (X) 

where the decision S*(X) provides the best return for that 
particular input value X. If X were to change the values of 
S*(X) and the return function should also change. 

Now consider a multistage problem as shown in figure 
2.2. This is the standard recursive DP structure with back- 
ward numbering where stage N is on the left and the stage 1 
on the right. The optimal return for stages i,i-1,...,1 is 
represented by Fi(Xi). The solution technique solves 
succeeding stages from right to left providing the optimal 
return FN (XN) up to the current stage under consideration. 
In the single stage problem Y was the output state variable, 
but now Y serves as the input to the next stage so 
Xi- 1 =Y=Ti (Xi,Si) . Thus, in an N stage problem, the output 
from stage i becomes the input to stage i- 1 for i=1,...,N. 

In figure 2.2 

Xi is the state variable 
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I 




Si 



SI 




Xi- 1 




XO 






Figure 2.2 Example of an N Stage Decision Problem. 

Si is the decision variable 
Ri is the stage return function 
Ti is the stage transformation function 

* after an Xi or Si variable indicates the optimal 
value for that stage, 

where i = 1,...,N is .the stage index. For each stage a deci- 
sion Si* * (Xi) will be made that yields the optimal return 
Fi(Xi) as a function of the present stage return 
R (Xi,Si* (Xi) ) and the optimal return Fi-l(Xi-l) for the 
stages 1 through i-1. The recursive equation for a mini- 
mization procedure would be, 

Fi (Xi) =min (Ri (Xi, Si) +Fi- 1 (Xi- 1) ) i = 1,...,N, (2.7) 

Si 

wher e 



Xi- 1=Ti (Xi,Si) =Xi-CiSi i = 1,...,N. (2.3) 

The stage transformation is simply the budget minus the 
guantity purchased multiplied by its cost. This recursion 
assumes the overall return function is the weighted sum of 
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the individual return functions for all replacement parts. 
The equation for this relationship is, 

MSRT=£Ei (Xi,Si) i=1,...,N. (2.9) 

The problem is to make a series of decisions SN...S1 that 
will optimize the overall return function while remaining 
within the budget constraints. 

For the previously discussed provisioning problem, using 
the different replacement parts as the stages, the state 
variable Xi is the resource, i.e., the dollars remaining at 
each stage. The decision variable Si is the number of parts 
purchased in stage i. The stage return function is the MSRT 
which must be computed from four operating characteristics 
for each part, (essentiality, procurement lead time, demand 
rate, performance as a function of quantity procured). The 
following equation represents the return function at each 
stage , 

Ri (Xi, Si) =Ei*Li*Ti*Zi (Si) /_Ei*Li*Ti (2. 10) 

i = 1 , . . . , N. 

The Standard DP approach computes the function Fi(Xi) 
from equation 2.7 recursively for i=1,...,N and for all 
possible values of Xi. At each stage i this is done by 
examining all feasible decisions for each Xi, computing the 
return for each decision, and then adding the corresponding 
optimal return Fi(Xi-l) from the remaining stages to arrive 
at the optimal return for the present Xi. Then for all Xi 
values in the stage, the best total return is chosen and 
assigned to Fi (Xi) . At the Nth stage the process is 
complete and FN(XN) represents the optimal return for the 
problem. This process is efficient for small problems, 
however it becomes computationally more difficult as the 
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number of stages and the number of possible state variable 
values increase in each stage. 

For an example of the computational burden, consider a 
very simple problem with a budget 3, a constant increment K 
for the state variable Xi within each stage, and a max of D 
decisions for each state variable value. In this case the 
number of evaluations for one stage is equal to (B/K) *D. 
Assuming the transformation function and the recursive equa- 
tion for the optimal return equation require only one addi- 
tion each, then there are 2*(B/K)*D additions per stage. 
Each state variable value will require D- 1 comparisons, 
hence, there are (B/K)*(D-1) comparisons per stage. The 
total number of operations (additions plus comparisons) 
required for all stages may be represented by, 

N ( (2*B*D) /K) + (B* (D-1) /K) , (2.11) 

or simplified this becomes 

(N *B/K) * (3D- 1 ) . (2.12) 

This example is actually a gross over-simplification of the 
computational problem. It is possible the return function 
could be very complex requiring many operations for each 
decision at each Xi value through the range of the budget. 
However, for a numerical comparison, consider a 25-stage 
problem with a budget of 74000, a state variable increment 
of 1000 and a maximum of 100 decisions. This leads to 
553150 operations as shown in equation 2.13, 

( (25*7400 0) /I 0 00) * ( (3* 100 ) - 1 ) =553 , 1 50 . (2. 13) 

This is the required number of additions and comparisons for 
a 25 stage problem. For today's high speed computer an 
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addition operation at machine language level takes about 
twenty-four machine cycles while a comparison takes 
eighteen. Using twenty cycles as an average and a machine 
speed of two million cycles per second, the above number of 
comparisons and additions would take approximately nine 
seconds. This time estimate is for an incredibly simple 
problem that does not take into account any of the adminis- 
trative overhead, objective or constraint computations or 
increased calculations required by complex functions. By 
setting K=10, which is equivalent to reducing the stage 
increment from 1000 to 10, equation 2.14 shows that the 
number of operations will increase to 

( (25*7400 0) /I 0)* ( (3*100) - 1) =55,315, 00 0. (2. 14) 

Hence, it is easy to see the limitation of standard dynamic 
programming. While it gives integer solutions, the computa- 
tional burden is considerable even in simple problems. 

C. A DYNAHIC PROGRAHHING VARIATION 

The computational workload of the full dynamic program- 
ming procedure makes it undesirable for large problems. A 
technique that will be referred to as •FUNNELING 1 because of 
its control of the resource boundaries between the stages, 
was developed and applied to a standard recursive DP 
computer program. The result is a shortened procedure that 
provides a locally optimal integer solution to problems that 
can be formulated in the DP construct. 

For any optimal solution to an N-stage DP problem there 
is a series of decisions Si* for i=1,...,N. Corresponding 
to this solution there is a sequence of Xi* values for 
i=1,...,N. These optimal state variable values Xi*, repre- 
sent the best quantities of resource available for 
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consumption at that stage and are not known in advance. If 
it were possible to find an approximate value for Xi* in 
each stage, the optimal solution then could be localized 
with an iterative process. That is, rather than compute 
F(X) for all X values there could be, a ’window’ of 
interest, centered around the approximate Xi* value in each 
stage. The computations would be performed on the limited 
range of Xi values in that window thereby reducing the 
required workload significantly. This concept of control- 
ling computational boundaries around an approximate optimal 
Xi* value is the foundation for the D? variation. 

The problem initially is how to localize the regions of 
the best Xi. This is solved by the first run of the funnel 
program using a large state variable increment. The 
resulting Xi values are then used as an approximation to the 
initial Xi* values i=1,...,N. A subroutine of the funnel 
program then constructs a window around the area of interest 
for further refinement. The first estimate of the Xi* 
values is quite rough and the decisions, especially the Si 
which have costs less than the initializing increment, may 
change significantly as the Xi are examined with a lesser 
state variable increment. 

The second program run is the first refinement of the 
initial solution. Using the large increments the program 
will have chosen the best return available though it may not 
be close to the optimal value. The increased sensitivity of 
the reduced state variable increment will allow considera- 
tion of more Xi values in the window permitting the opti- 
mization process to move closer to the best return. It is 
conceivable that any Xi may change as much as the window of 
its stage will allow. This much variation in the decisions 
indicates that the best solution lies outside the window and 
the boundaries must be changed accordingly. This situation 
may be present in several stages for any program iteration. 
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Due to the characteristics of the DP procedure, a change in 
the state variable consumed at stage N will cause the 
quantity of resource available at stage N-1 and succeeding 
stages to change. The concept is that repetitive runs of 
the funneling program will allow the computational bounds 
for each stage to localize and center upon their respective 
Xi*. When each of the stages has chosen its Xi*, the Si* 
for i=1,...,N will represent the best solution for a given 
objective and budget. This solution will be a local and 
possibly a global optimum. Further discussion of the 
program specifics may be found in Chap 3. 

Additional computations may be saved in the funneling 
program by ranking the stages, largest to smallest, by the 
cost of their respective parts. Purchasing the most expen- 
sive items first will reduce the magnitude of the budget and 
simplify the computations in the remaining stages. The 
amount of computation depends on three things, 1) the width 
of the window for each stage, 2) the number of state vari- 
able increments within the window and 3) the number of deci- 
sions available for each state variable increment. 

- Consider again the computational example of section 3, 
where the number of operations is estimated for the dynamic 
programming solution procedure. To use that analysis for 
the funnel procedure, B in equation 2.12 is replaced by the 
width of the window W. This width will vary between stages 
but this example will assume Wi is constant and equal to 
10000. Setting the remaining factors the same, K = 10, S=100, 
and N=25, the number of operations required from equation 
2. 1 2 is, 

{ (25*10 000) /I 0) * ( (3*100) -1) =7,475,00 0. (2. 15) 

This is roughly 13.57c of the full DP requirement computed in 
equation 2.11, a marked improvement in efficiency. There 
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are other ways to make the algorithm more efficient and save 
computations. Further discussion of the details of the 
modified DP procedure may be found in chapter 3. 
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III. THE FUNNEL PROGRAM 



This chapter discusses the operation of the standard DP 
program, the modifications required to employ the funneling 
concept and some of the implementation problems encountered 
in the process. Specific subroutine operations are 
explained in the second section. 



A. MODIFICATION FOR THE FUNNEL PROCEDURE 



The DP variation, called the Funnel program, is a modi- 
fied version of a standard DP computer program named DP5. 
The modified program will be referred to as DP6 for the sake 
of clarity in the following discussion. The DP5 program as 
described in [Ref. 3] has four subroutines, ST3P.ET, STORE, 
TRANFM and DLIMIT. The following list describes what func- 
tion each subroutine performs in the DP process. 



1) STORE 



- enters the constants to be used in the other 
subroutines and the main program. 



2) DLIMIT 



defines the range of the decision values Si 
for any value of Xi in stage i. 



3) STGRET - defines and computes the return function for 

any values of Xi and Si in stage i, Ri(Xi,Si). 

4) TRANFM - defines and computes the transformation 

function for any values of Xi and Si in stage 
i, Ti ( Xi, Si) . 



Providing DP5 with the specific input data as required in 
[Ref. 3], the program will compute the optimal return by the 
standard recursive DP procedure. The following figure shows 
the flow diagram for the general solution technique. 
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Figure 3. 1 
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Figure 3.2 Flow diagram of the DP6 procedure. 
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When there is no previous solution, the beginning data file 
may be constructed in the correct format as defined by the 
user instructions located in the comment section of DP6, see 
appendix. This input file, as described, may be used 
directly with the program. The next step in the figure is 
the first iteration of the Funnel process utilizing the 
DLIMIT, TRANFH, and STGRET subroutines. When the calcula- 
tions are complete the results must be examined to see if 
the process has found the optimal solution, i.e., 'is it 
DONE'. The requirements for the 'DONE' condition are 1) an 
all integer solution, 2) no window boundaries violated and 
3) the budget remaining must be less than the smallest cost 
for items to be purchased. If all of these conditions are 
not met the ADJOST subroutine is activated and the input 
data file is manipulated to correct the problem areas. With 
each program iteration this corrective procedure is repeated 
until all the conditions are met and an optimal solution is 
produced. An example of the Funnel program output is shown 
is figure 3.3. 



THE OBJECTIVE IS MSRT 
CONSTHAINT ON COST 

THE DESCRIPTION WHICH APPEARS BEICW FOR STAGE 1 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE 
THE PROBLEM IS TO MINIMIZE A 3 STAGE PROCESSXN 
IS TO BE CHOOSEN OPTIMALLY BETWEEN XN=2- 05000C+02 
AND XN=2- 0500D + 02 

OPTIMAL IN = 2.05000D«-02 OPTIMAL BETO RN=5. 30 206 D+ 0 0 
N XN SN XN-1 

3 2-0500D+02 1.00D+01 5.5000D+01 

2 5.5000D+01 2. 00D+Q0 3.50C0E»01 

1 3- 5000D*0 1 7.00D+00 0-0 



Figure 3.3 Example of DP6 Output. 
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The objective and constraint are stated in the first two 
lines. The description that follows gives the number of 
stages, the budget quantity, the optimal return and detailed 
stage results. Specifically delineated for each of the 
stages are the input budget quantities, the decisions and 
the output budget quantities. For this particular example 
the budget remaining at the- end of the problem is zero. 

Two ancillary functions were added, one to the main 
program and one to a subroutine to provide information to 
the user. The first is a timing function that provides the 
program run time. The second is a counting function that 
presents to the screen and stores in a file the numbers of 
the stages in which the optimal -return values are too close 
to the computational boundaries. 3oth functions have 
provided valuable insights that will be discussed in Chapter 
5. 

B. THE SUBROUTINE ADJUST 

The subroutine ADJUST is the heart of the funneling tech- 
nique. To gain the computation reduction desired, this 
subroutine manipulates the windows around the ’current* 
state variable values, Xi' , in each stage for each program 
iteration. Initially the upper and lower limits HXi’ and 
Lxi' on the window at stage i are defined by, 

HXi* = Xi ' +K*Ci+ 1 , (3.1) 

and 

LXi’ = Xi’-K*Ci+ 1 , (3.2) 

where K is a constant. These bounds are illustrated in fig 
3.4. Further adjustment of these limits may be necessary to 
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assure that one of the Xi values considered in the computa- 
tion of F (Xi) is Xi*. At stage i in the computation of 
F(Xi) the optimization considers first Xi=LXi, then Xi is 
incremented by Ci to give the next value, etc. Thus if LXi 
is not properly set, the value Xi=Xi' will never be 
considered. 



« HXi 



K x Ci+1 



Xi 




LXi 



LXi->- 



Figure 3.4 Example of a Hindow on a Stage i 
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This adjustment of LXi' is made by decrementing from the 
temporary stage value, Xi' , by the cost Ci for that stage. 
When this decrementing process first yields a value less 
than or equal to LXi’ that value is set as the new lower 
limit on the window and is called LXi. A similar adjustment 
is made to the upper limit by incrementing from Xi' until 
the value HXi' is equalled or exceeded. The new upper limit 
is called HXi. This process is repeated for al‘l the stages 

until the funnel is formed. In figure 3.5, the funnel may 

be seen for N stages of a provisioning problem. The budget 
quantity B is represented by the vertical lines for all 
stages. The upper and lower bounds are represented by HXi, 
LXi for i=1,...,N where i is the stage index. In each stage 
the bounds are at least a distance equal to the product of 
the constant K and Ci+1 from the Xi’ value. Stages where 
the bounds are further from the center indicate the addi- 
tional adjustment for compatability between the successive 
Xi’ values. 
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Figure 3.5 A Funnel Diagram for N Stages of a Problem. 



The reason for multiplying the succeeding stage costs by 
the constant K is to allow the DP procedure some computa- 
tional freedom in each of the stages while searching for the 
Xi and Si*(Xi) values. During the computational process, 
should a stage choose an Xi • that is within Ci+1 of the 
computational bound, the boundary is said to be violated. 
This is because the DP procedure may actually prefer a value 
outside the bound but is unable to reach it. Figures 3.6 
through figure 3.9 show an example of the stage boundary 
manipulation as the iterative process of the DP6 program 
searches for the optimal solution. 



I 

i 

\ 

1 



t 



i 

—200 =HXi 



i-150 -Xi ’ 



—100 =LXi 



i-1 



-pl70 =HXi- 1 

' -^L2Q.=Xi-l ' 
jl70 =LXi- 1 



i-2 



T 1 20 =HX i- 2 

j; 70 =Xi-2 ' 

\ 

X 20 =LX i - 2 



l 



Figure 3.6 The Initialization Run. 



The initializing run depends on the starting condition for 
the problem. If the ENTER function is used, the state vari- 
able increments will represent the costs of the respective 
stages. Figure 3.6 shows a starting condition of no prior 
knowledge; hence, the state variable increment is a 
constant, in this example 10, for all stages and of egual 
number on either side of the current Xi' estimate. 

A second run of the program provides an opportunity for the 
DP process to evaluate the estimated optimal returns with a 
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i- 1 



_ 1 70=HXi- 1 



120=HXi-2 




i 1 0 0 = LXi 



Ci= 10 



f =Xi- 1 ' 
-j_7Q=LXi = l _ 



Ci- 1 = 5 



5 0 = X i - 2 ' 

i: 

2 0 = LXi- 2 



Figure 3.7 The Second Run. 



smaller state variable increment or a revised boundary value 
and to adjust the window around a new Xi • accordingly. 
Figure 3.7 shows the funnel after the second run tut before 
the windows are adjusted. The program has selected Xi*=140, 
Xi-1’=70, and Xi-2’=50. These are the current estimates of 
the optimal Xi ' values at each stage. Notice that stage i- 1 
has a boundary that is violated, i.e., the Xi-1' for the 
stage is within Ci=10 of the boundary. In this case the 
optimization function may prefer an Xi-1 outside the window 
but the boundary prevents its selection. Also in the second 
run the state variable increments have been changed to the 
cost of the item associated with the stage. In most stages 
these values are smaller then the initial increment making 
the return functions more sensitive. This increase in 
sensitivity should improve the optimization selection. 

The violation of the lower boundary in stage i- 1 causes 
the iteration to fail the conditions of the DONE checkpoint 
of figure 3.2. Consequently, the ADJUST subroutine is acti- 
vated again to manipulate the boundaries producing the 
results shown in figure 3.8. 
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Figure 3.8 The Window Adjustments. 
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Figure 3.9 The Third Run. 
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Figure 3.9 shows the funnel after the third run. The 
program has selected Xi’ = 140, Xi-1'=65 and Xi-2'=46 for the 
current estimates of the Xi 1 values at each stage. Notice 
that the Xi-1* value is just outside the previous LXi-1, 
which was 7 0. In this example the Funnel program manipu- 

lated the boundaries of the windows to allow the DP process 
to choose different Xi' values for the three stages, i, i-1, 
i-2. Subsequent runs of the program will determine the Xi* 
for each stage. The process is 'DONE 1 when the remaining 
budget <= minimum Ci, an optimal integer solution is 
provided, and there are no boundaries violated. 

C. IMPLEMENTATION PBCBLEHS 

From the characteristics of the funneling process, there 
are problems that arise when transitioning between stages. 
When computing the maximum and minimum value of the decision 
to be made for each Xi in a stage i, the upper and lower 
bounds and item cost of stage i- 1 must be considered. 

In figure 3.10, the left side is the window for stage i, the 
right side is the window for stage i-1. For the stages i, 
i-1, the highest values of the state variables are repre- 
sented by HXi and HXi-1 while the lowest values are repre- 
sented by LXi, I Xi-1 . In the program DP6, the two variables 
DLO X and DHIGH represent the minimum and maximum decisions 
available for the state variable in each stage. It is 
essential that DHIGH for any Xi value in the window of the 
ith stage be of such a magnitude that it will not force the 
Xi-1 value below the lower bound LXi-1. Additionally, DLOW 
must be of such magnitude that for any Xi in the window of 
the ith stage, it will not force the Xi-1 value to violate 
the upper bound HXi-1. If either of these decision quanti- 
ties is incorrect at stage i the resource quantity will be 
outside the window at stage i-1. The derivation of the 
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Figure 3.10 Boundary Manipulation. 

formulas for these values follows. The transformation func- 
tion as represented by equation 2.8 of Chap 2 is, 

Xi-1=Xi-CiSi i = 1,...,N. (3.3) 

If this transformation is to take place from a window in 
stage i to a window in stage i-1 the following conditions 
must be true. 

LXi-1 <= Xi-Ci*SHIGH, (3.4) 

where SHIGH is the maximum number of items to be purchased 
and 



Xi-Ci*SL0W <= HXi-1, (3.5) 

where SLOW is the minimum number of parts to be purchased. 
By algebraic manipulation a relationship for SLOW and SHIGH 
may be established. 
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SHIGK <= (Xi-LXi- 1 )/Ci, 



(3.6) 



and 

SLOT? >= (Xi-HXi-1) /Ci. (3.7) 

The expressions SHIGH AND SLOW are equivalent to the DHIGH, 
DLOW required in DP6 and are coded in the subroutine D LIMIT 
to control the boundary calculations. It is conceivable 
that at the extreme state variable values of Xi the computed 
boundaries could become negative or larger than the assigned 
budget. Thus, the decisons are further limited by DHIGH= 
min (SHIGH, 2 00) and DICW= nax(0, SLOW) , where 200 is an arbi- 
trarily chosen value. 

The initialization for the problem is either from input 
decisions or else a large state variable increment is used. 
After the subroutine ADJUST manipulates the windows the 
first time, the state variable increments are decreased or 
the boundary conditions are changed to allow examination of 
the new possible Xi values for the best return. The subrou- 
tine DLIMIT automatically adjusts the DHIGH and DLOW values 
to compensate for the changes in the state variable incre- 
ment. With this closer scrutiny, the optimal return value 
may change. The ADJUST subroutine moves the boundaries up 
or down accordingly when an optimal value is within one 
increment of the bounds. Recall that if the optimal return 
value is within one increment of window bound, it may be 
that the DP process would prefer a value outside the window. 
Hence, the boundary is recorded as violated, and the window 
is moved for the next iteration. This process is repeated 
until all the optimal return values are centered within the 
windows for all stages. 

To save computations and avoid computing Xi values that 
are not compatible between stages, the state variable 



4 1 



integer, 
cons idered. 



increment for each stage is the cost of the item to be 
purchased in the previous stage. A problem arises in that 
the original DP5 program may provide solutions that are not 
This is possible because in the stage being 
the best return may require an Xi ' value from 
succeeding stages that has not been computed. DP5 handles 
this by interpolating linearly between the state variable 
values it does have to arrive at an estimate for the one 
required. In DP6 tc avoid this approximation and obtain 
integer solutions, the subroutine ADJUST, (as discussed in 
section A), manipulates the lower and upper window values to 
ensure the previous stage optimal return values are avail- 
able for the Xi being considered. This adjustment process 
is also repetitive and is done in conjunction with the 
boundary control. 

Another problem arises from the effort to save computer 

If there is a wide range of costs, the 
the state variable increment within 
Initially the ADJUST subrou- 
tine creates the state variable windows for the individual 
stages using a constant multiplier that may be set by the 
user. To remedy the problem of the extreme range of the 
-costs of items for the provisioning problem being consid- 
ered, ($10000.00 to .50), the multiplier was changed to a 
variable that depends directly on the item cost. This was 
done to decrease the program repetitions required to reach 
an optimal solution. Under the constant multiplier method, 
a change in the purchase quantity of a larger cost item, 
(i.e.$1500), could require many program repetitions to 
re-allocate the resource available for the lower cost items, 
(i. e. $10,5, 3.75, etc). The smaller window range, due tc the 
constant multiplier, hampers the re-allocation effort. The 
variable multiplier allows greater compensation for lower 
cost items to facilitate a more rapid change if required. 



computation time, 
use of the item cost as 
a stage becomes troublesome. 
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IV. ANALYSIS 



The primary topic of this chapter will be the comparison 
of the Marginal Analysis and Funneling solution techniques 
for the previously stated problem. The discussion will 
cover the results of both programs and the solution differ- 
ence. The standard dynamic program procedure will net be 
addressed because of the computational advantage of the DP 
variant and the similiarity of the two techniques. 

Tor simplicity and ease of comparison a small three 
stage problem will be defined and solved using both solution 
procedures. As in the original provisioning problem the 
objective is to minimize the MSRT with a constraint on cost. 
Three different budget quantities will be examined to demon- 
strate the differences in the solution techniques. The 
following table is a list of the input parameters which are 
cost Ci, demand rate Li, essentiality Ei, and procurement 
lead time Ti. For this particular example Ti is one year. 



TABLE XI 

Sample Problem Input Parameters 



I DEMAND COST 

1 5.000 5.00 

2 1.000 10.00 

3 10.000 15.00 



T ESSEN 

1.00 3.0 
1.00 2.0 
1.00 1.0 
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For the first budget quantity of 3195.00 the Marginal 
Analysis program returns the following output. 







TABLE III 

Marginal Analysis Program Return, Budget=$ 1 95. 00 

*********************************************** 

* METHOD: FIXED BUDGET = 195-00 * 

*********************************************** 

*********************************************** 



STEP 


ITEM 


SI 


BATIO 


BB 


1 


1 


1 


0-480808496 


190-00 


2 


1 


2 


0-365659833 


185- 00 


3 


1 


3 


0-260618091 


180-00 


4 


1 


4 


0- 172421157 


175-00 


5 


1 


5 


0- 105280340 


170-00 


6 


2 


1 


0.0735758543 


160-00 


7 


3 


1 


0.0600002967 


145-00 


8 


1 


6 


0-0591956675 


140.00 


9 


3 


2 


0-0533369593 


125-00 


10 


3 


3 


0-046688754 1 


1 10. 00 


1 1 


3 


4 


0.0400909968 


95.00 


12 


3 


5 


0-0336193480 


80-00 


13 


1 


7 


0.0306577347 


75.00 


14 


3 


6 


0-0273999237 


60-00 


15 


3 


7 


0.0216008648 


45- 00 


16 


2 


2 


0.0207276568 


35-00 


17 


3 


8 


0-0164023377 


20.00 


18 


1 


8 


0.0146531016 


15.00 


19 


3 


9 


0.0119544677 


0.0 



****************************************** 


MINIMIZE 


MSET 


EXCESS: 0-0 


I 


SI 


MS BT 


1 


8 


1-302 


2 


2 


10.396 


3 


9 


12. 176 




OVEEALL 


6-003 



REALIZED PERFORMANCE = 6.0033 

BOUNDS ( 6-0033, 6.0033) 



The decision results are, S1=8, S2=2, S3=9 and the return is 
Z= 6.0033, with no budget remaining. Viewing the purchase 
sequence in table III it may be seen that the Marginal 
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that 



Aiialysis program has reached a solution uninterrupted, 
is without implementing the heuristic developed in [Ref. 1], 
Consequently, for this budget quantity the solution is 

optimal. Notice another indicator of the optimal solution 
in the last line of table III. The bounds described are the 
possible upper and lower values for the return. When these 
quantities are the same value, the solution is the best 
available. 



TABLE IV 

Funnel Program Results, Budget=$ 195. 00 



THE OBJECTIVE IS HSBT 
CONSTfi AI5T ON COST 

THE DESCRIPTION WHICH APPEARS BEICW FOR STAGE 1 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE 
THE PROBLEM IS TO MINIMIZE A 3 STAGE PROCESSXN 
IS TO BE CHOOSER OPTIMALLY BETWEEN XN=1 - 9 5 00D *0 2 
AND XN= 1. 9500D + 02 

OPTIMAL XN=1- 95000D+02 OPTIMAL EETURN=6.00332D+00 
N XN SN XN-1 

3 1- 9500D*02 9.00D+00 6.0000D + 01 

2 6. 00 0 0 D* 0 1 2- Q0D+-QQ 4.0000D+01 

1 4. 00 00D+ 0 1 8.00D + 00 0-0 



The same budget and parameters entered in the Funnel 
program produce the results in table IV. The output is in 
the same format discussed in Chapter 3. The initial budget 
quantity, XN and the optimal return are shown in the line 
just above the individual stage results. The variable Xi is 
the budget quantity at the beginning of the stage. Si is the 
decision or quantity purchased, and Xi-1 is the budget at 
the end of the stage. Recall that item 3 is the most expen- 
sive and item 1 the least expensive. The optimal solution 
for the $195.00 budget is S1*=8, S2*=2, and S3*=9 with a Z 
value equal to 6.0033 and no budget remaining. This answer 
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matches exactly the Marginal Analysis program results shown 
in table III, confirming that the Marginal Analysis solution 
is optimal. 



TABLE V 

Marginal Analysis Program Return, Budget=$200 . 00 

*** * * * ** * * * * * ** ** *r** '* * **n* * '* ** F * * 4 '** ♦ * ****** *~* * * 

* METHOE: FIXED EODGET = 200.00 * 

*********************************************** 

******************************** *************** 



STEP 


ITEM 


SI 


RATIO 


BR 


1 


1 


1 


0.480808496 


195.00 


2 


1 


2 


0.365659833 


190.00 


3 


1 


3 


0.260618091 


185.00 


4 


1 


4 


0. 172421157 


180.00 


5 


1 


5 


0.105280340 


175.00 


6 


2 


1 


0.0735758543 


165.00 


7 


3 


1 


0.0600002967 


150.00 


8 


1 


6 


0.0591956675 


145. 00 


9 


3 


2 


0.0533369593 


130.00 


10 


3 


3 


0.0466887541 


1 15.00 


1 1 


3 


4 


0.0400909968 


100.00 


12 


3 


5 


0-0336 193480 


85. 00 


13 


1 


7 


0.0306577347 


80. 00 


14 


3 


6 


0.0273999237 


65.00 


15 


3 


7 


0.0216008648 


50. 00 


16 


2 


2 


0.0207276568 


40.00 


17 


3 


8 


0.0164023377 


25.00 


18 


1 


8 


0.0146531016 


20. 00 


19 


3 


9 


0.01 19544677 


5.00 


20 


3 


10 


0.0083406679 


-10.00 


20 


1 


9 


0.0064818673 


0.0 



*********************************************** 

MINIMIZE MSRT EXCESS: 0.0 



1 
1 

2 
3 



SI 

9 

2 

9 



MSRT 
0. 514 
10.396 
12. 176 



OVERALL 5.565 



REALIZED PERFORMANCE = 5.565 

BOUNDS ( 6.0033, 4.3120) 



For the second example the budget quantity 
was used. The Marginal Analysis program returns 



of $200.00 
the output 
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in table V. The results are S3=9, S2=2, and S1=9 with a Z 
value of 5.5652 and no budget left over. Examining the 
detailed output of the Marginal Analysis program, table V, 
it may be seen that the purchase of the twentieth item 
drives the budget negative. The program implements the 
heuristic, re-evaluating the benefit to cost ratios, then 
chooses the number one item for the next purchase. Notice 
that in the bottom line of table V the upper and lower 
bounds are not equal. This is a good indication the solu- 
tion may not be optimal and should be verified. 



TABLE VI 

Funnel Program Results, Budget=$200. 00 



THE OBJECTIVE IS tISRI 
CONSTRAINT ON COST 

THE DESCRIPTION HHICR APPEARS EEICW FOR STAGE 1 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE 
THE PROBLEM IS TO MINIMIZE A 3 STAGE PROCESSXN 
IS TO BE CHOOSEN OPTIMALLY EETHEEN XN=2- 00C0D+02 
AND XN=2- 0000D + 02 „ 

OPTIMAL XN=2- 000JD+02 OPTIMAL R£TUBN=5. 565 19D+00 
N XN SN XN-1 

3 2. OOOOD+02 9. QOD«-00 6.500CD + 01 

2 6- 5000D + 0 1 2.0CD + CQ 4-5000D+01 

1 4.5000D+01 9-OOD+OO 0-0 



The same budget and parameters entered in the Funnel 
program produce the results given in table VI. In this case 
the answer is again identical, S1=9, S2=2, and S3=9 with the 
optimal return' Z= 5.56519. Incidentally, the heuristic in 
the Marginal Analysis program makes the correct choice at 
the twentieth step when it selects the item with the second 
highest benefit to ccst ratio. 

For a third case, a budget of $205.00 is entered in both 
programs. The results for the Marginal Analysis code are 
shown in table VII. 
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TABLE VII 



Marginal Analysis Program Results, Budget=$205. 00 

*********************************************** 

* METHOD: FIXED EUDGET = 205-00 * 

*********************************************** 

*********************************************** 



STEP 


ITEM 


SI 


RATIO 


BB 


1 


1 


1 


0-480808496 


200.00 


2 


1 


2 


0.365659833 


195.00 


3 


1 


3 


0. 260618091 


190.00 


4 


1 


4 


0. 172421157 


185.00 


5 


1 


5 


0. 105280340 


180-00 


6 


2 


1 


0.0735758543 


170.00 


7 


3 


1 


0.0600002967 


155.00 


3 


1 


6 


0.0591956675 


150.00 


9 


3 


2 


0.0533369593 


135.00 


30 


3 


3 


0.046683754 1 


120.00 


11 


3 


4 


0.04C0909968 


105. 00 


12 


3 


5 


0.0336193480 


90.00 


13 


1 


7 


0.0306577347 


85.00 


14 


3 


6 


0.0273999237 


70.00 


15 


3 


7 


0.0216008648 


55. 00 


16 


2 


2 


0-0207276568 


45.00 


17 


3 


8 


0.0164023377 


30- 00 


18 


1 


8 


0.0146531016 


25.00 


19 


3 


9 


0.0119544677 


10.00 


20 


3 


10 


0.0083406679 


-5.00 


20 


1 


9 


0.0064818673 


5.00 


21 


1 


10 


0.0026624922 


0 . 0 



*********************************************** 

MINIMIZE MSBT EXCESS: 0.0 



1 
1 

2 
3 



SI 

10 

2 

9 



MSB! 
0. 190 
10.396 
12. 176 



OVERALL 5.385 



REALIZED PERFORMANCE = 5.3852 

BOUNDS { 6.0033, <1.3120} 



The solution is S1 = 10, S2=2 , and S3=9 with a return value of 
Z= 5.3852. The upper and lower bounds are not equal indi- 
cating the heuristic has been applied and the solution 
should be verified. Examining the output from the Marginal 
Analysis program, table VII, it may be seen the budget was 
again driven negative at the twentieth item. The heuristic 
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was activated purchasing two of the number one items to 
drive the budget to zero. 

Entering the same input and budget parameters in the 
Funnel program produces the results shown in table VIII. 



TABLE VIII 

Funnel Program Results# Budget=$205. 00 



THE OBJECTIVE IS MSRT 
CONSTRAINT ON COST 

THE DESCRIPTION WHICH APPEARS BELOW FOB STAGE 1 
APPLIES TO STAGE 1 THRU STAGE 3 INCLUSIVE 
THE PROBLEM IS TO MINIMIZE A 3 STAGE PROCESSIN 
IS TO BE CHOOSEN OPTIMALLY BETWEEN X N=2- 05000 D ♦ 02 
AND XN=2. 05Q0D+02 

OPTIMAL IN=2«050Q0D*02 OPTIMAL BETURN=5. 30246D+00 
N XN SN XN-1 

3 2- 0500D+02 1.00D+01 5-5000D+01 

2 5. 5Q00D*0 1 2. OOD *• 0 0 3.50C0D*01 

1 3- 5000D*- 0 1 7.00D+00 0.0 



The answers are not the same. The Funnel program results 
indicate the best part selection is S1 = 7, S2=2, and S3=10 

for an optimal return of Z= 5.30246. This selection 
provides an MSET that is roughly .08 better than the 
marginal analysis method. 

Recall from section A of Chapter 2 that the Marginal 
Analysis program has a function called KIT that will deter- 
mine the return from a given decision set. Entering the 
optimal solution from the Funnel program into the KIT func- 
tion produces the result in table IX. The resulting objec- 
tive value is identical to the value computed by the DP 
variant. This shows that the computational procedures are 
the same for both programs. 

Now consider the definition of a much larger problem. 
It is of the same type as the previous examples, an 
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TABLE 


IX 


Kit Function Results, Budget=$205. 00 


EVALUATION 


OF SPECIFIED ALLOCATION..-. 


I 


SI 


MSRT 


t 


7 


3.085 


2 


2 


10.396 


3 


10 


7.610 




OVERALL 


5-302 


**** **************************** *************** 


' 




1 



objective of minimizing MSRT with a constraint on cost. The 
input parameters cost Ci, demand rate Li, essentiality Ei, 
and procurement lead time Ti, are shown in in table X. The 
budget is $74825.00 with twenty five stages and costs 
ranging from $10000.00 to .50. Due to the length of the 
detailed results an abbreviated output of the Marginal 
Analysis solution is shown in table XI. The return is 2. 853 
with this decision set. Note in the last line of the table, 
the split in the bounds indicates the heuristic was imple- 
mented and the solution may not be optimal. Examining the 
abbreviated output of the Marginal Analysis program. Table 
3.4 shows that the budget was driven negative ten iterations 
from the end. Entering the same parameters and budget quan- 
tity for bcth starting conditions of the Funnel program 
produces identical results shown in table XII. Comparing 
the results in tables XI and XII it may be seen that the 
return and decision set for the problem are the same for 
both programs. Interestingly, the Marginal Analysis program 
was able to correctly select ten items after the heuristic 
implementation to match the optimal solution produced by the 
Funnel program. The heuristic employed by Richards and 
McMasters appears to have been a good choice, although the 
solution must still be verified. 
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TABLE X 

Input Parameters for a Larger Problem 






I 


DE9AND 


COST 


T 


ESSEN 


1 


2-000 


0.50 


1.90 


1-0 


2 


17.000 


3-00 


2-20 


3-0 


3 


0-250 


3.75 


2-50 


3-0 


4 


50-000 


5.00 


2.60 


1.0 


5 


1.000 


10.00 


1-60 


3.0 


6 


1- 500 


12-00 


2.50 


1.0 


7 


18-000 


15.00 


1-50 


3-0 


8 


3-000 


20.00 


2-50 


1.0 


9 


0.250 


25.00 


1.50 


1.0 


10 


1.000 


25.00 


2.00 


1-0 


1 1 


2-000 


25.00 


1-50 


2.0 


12 


0.330 


50.00 


1-50 . 


3-0 


13 


0-100 


50-00 


2-0 0 


3-0 


14 


1.000 


61.00 


2-00 


2-0 


15 


7.000 


80-00 


2.30 


2.0 


16 


1-000 


91.00 


3-70 


3-0 


17 


2. 500 


150.00 


1-75 


3-0 


18 


20.000 


210.00 


1.70 


3.0 


19 


3.000 


450.00 


2.70 


2-0 


20 


0.500 


500.00 


1.75 


1.0 


21 


1.500 


710.00 


3. 00 


2.0 


22 


6.000 


1000.00 


1- 80 


1.0 


23 


2-000 


1500.00 


2.50 


3.0 


24 


12-000 


1500.00 


2.00 


2-0 


25 


0-500 


10000.00 


1- 00 


1.0 
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TABLE XI 



i 



Marginal Analysis Program Results, Budget=$74825. 00 

* 1 ***"* ****** *********** * *** ****** * * *■** ********** 

* METHOD: FIXED BUDGET = 74825.00 * 

*********************************************** 

******************************** 4 ************** 



STEP ITEM SI RATIO EP 



1 


1 


1 


2-82236671 


74824.50 


2 


2 


1 


2. 141 17622 


74821.50 


3 


2 


2 


2.08235264 


74818-50 


4 


2 


3 


2.02352905 


74815.50 


5 

m 


2 

<• 

• 


4 

• 

• 


1-96470547 

• 

m 


74812.50 

- «• 

m 


440 


24 


21 


ol 0004 187953 


23861 50 


441 


2 


55 


0.0004181664 


2383-50 


442 


22 


9 


0.0003934940 


1383.50 


443 


19 


11 


0.0003861645 


933. 50 


444 


18 


40 


0-0003540027 


723-50 


445 


4 


154 


0.0003520278 


718.50 


446 


24 


22 


0.0003425654 


-781.50 


446 


12 


3 


0-0003396191 


668.50 


447 


15 


23 


0.0003359483 


588.50 


448 


7 


39 


0.0003278039 


573.50 


449 


8 


14 


0-0003014780 


553-50 


450 


20 


2 


0.0002939128 


53.50 


451 


4 


155 


0.0002806289 


48-50 


452 


2 


56 


0.0002605470 


45.50 


453 


10 


6 


0-0002369717 


20. 50 


454 


1 


12 


0.0002271836 


20-00 


455 


4 


156 


0.0002225503 


15-00 


456 


7 


40 


0.0002018004 


0-0 



REALIZED PERFORMANCE = 2.8532 

BOUNDS i 2.9796, 2-6762) 
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TABLE III 



Funnel Program Results, Budget=$74825. 00 

THE OBJECTIVE IS HSBT 
CONSTRAINT ON COST 

THE PROBLEM IS TO MINIMIZE A 25 STAGE PROCESS 
IN IS TO BE CHOOSEN OPTIMALLY BETWEEN XN=7. 48250D+04 
AND XN=7.48250D+04 



OPTIMAL 


XN = 7. 48250D+04 


OPTIMAL BETOBN 


=2.85315D*00 


N 


XN 


SN 


XN-1 


25 


7. 4825D+Q4 


0.0 


7. 4825D+04 


24 


7.4825D+04 


2. 10D+01 


4.33250+04 


23 


4.3325D + 04 


6. 00D*C0 


3. 4325D + 0 4 


22 


3. 4325D+04 


9. 00D*00 


2.53250+04 


21 


2- 5325D+04 


6- 00D* 0 0 


2. 1065D*04 


20 


2. 1065D+04 


2. 00D+C0 


2.00650+04 


19 


2. 0065D+04 


1. 10D + 01 


1.5115D + 04 


18 


1- 51 15D+04 


4. 00D+01 


6.7 150D + 03 


17 


6.71500*03 


8. 00D+00 


5.51500+03 


16 


5. 5150D+03 


8. 00D*00 


4.7870D+03 


15 


4.7370D+03 


2.30D+01 


2. 947CD*03 


14 


2. 9470D+03 


5. 000*00 


2-64200+03 


13 


2. 6420D+03 


2. OOD*0O 


2. 54200+03 


12 


2. 5420D+03 


3. 00D* CO 


2. 3920D+03 


1 1 


2.3920D+03 


7. 00D*Q0 


2. 2 170D + 03 


10 


2.2170D+C3 


6. 00D+C0 


2. 0670D + 0 3 


9 


2. 0670D+03 


2. OOD+OO 


2.0170D+03 


8 


2. 0170D+03 


1. 40D+01 


1.7370D+03 


7 


1.7370D+03 


4- 00D*0 1 


1. 13700*03 


6 


1. 1J70D+03 


9. 00D * 00 


1.02900+03 


5 


1. 0290D+03 


6. 00D+00 


9.6900D+02 * 


4 


9- 6900D+02 


1.56D + 02 


1.89000+02 


3 


1.8900D*02 


4. 00D+00 


1.7400D+02 


2 


1.74000*02 


5.60D+01 


6. OOOOD+OO 


1 


6. 0000D*00 


1.20D+01 


0.0 
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V. CONCLUSIONS 



This chapter will discuss the advantages and disadvan- 
tages of the two primary solution technigues examined in 
this paper. Various starting modes of the Funnel program 
will be looked at in an effort to develop an efficient 
approach to this type of provisioning problem. 

The Marginal analysis program provides a quick, simple 
and effective approach for generating a solution to the 
previously defined provisioning problem. However, the solu- 
tion is not guaranteed to be optimal because the program may 
implement the Richards heuristic to consume all available 
resources. The Funnel program provides the optimal solution 
to the same problem though it takes more computer time to 
operate. The DP variant has the capability to start the 
problem from two different conditions, 1) no prior knowledge 
and 2) prior knowledge or an estimate of the solution. For 
Case 2, utilizing the ENTER function of the Funnel program 
allows the rapid construction of the input data file used by 
DP6 to examine the Xi* values and manipulate the window 
boundaries as necessary. If the entered solution set is not 
optimal, the Funnel program will drive the decisions to the 
optimal decision set in subsequent iterations. Using the 
timing function mentioned in Chap. 3 a rough approximation 
of run time is available for each iteration of the program. 
For Case 1, no prior knowledge, the total estimated time for 
the $74825.00 problem is nearly 3000 seconds. For Case 2, 
prior knowledge, a near optimal starting solution for the 
same problem is driven to the optimal solution in approxi- 
mately 300 seconds. This 90% savings on computation time 
makes it clear that an approximation method as an input to 
the Funnel program would be a better approach to this type 
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of problem. It should be noted that the run time for the 
case 2 start is related to the accuracy of the entered deci- 
sion set. In this instance the entered decision set was 

specifically constructed to force the Funnel program to 
adjust the window boundaries in all stages at least once. A 
better solution than the one entered could produce a much 
shorter run time. 

There are tradeoffs in the Funnel Program for computa- 
tional efficiency. Variations of the multiplier that 
creates the window can drive the individual program run time 
directly. A window that is too large wastes computation 
time on resource values that are not needed. If the window 
width is too small, multiple iterations of the program may 
be required to compensate for the increased number of 
boundary violations in the effort to reallocate the resource 
consumed at the various stages. For the verification of a 
decision set provided by the Marginal Analysis program the 
multiplier could be set to 1.0 providing a very narrow 
window of resource to be examined. Due to the narrow width 



of the window, if the Xi* 
boundary will be recorded 



value changes in any stage the 
and displayed to the user as 
violated. This notification would indicate that the entered 
solution is not optimal for the problem. 

The computer time expended for the verification alterna- 
tive would be quite low even for large budget problems. In 
the 25-stage problem, for example, approximately one second 
is required to verify the correctness of the optimal solu- 
tion if it is entered initially. If the verification alter- 
native indicates the decision is not optimal, the time 
required to drive the entered solution to the optimal solu- 
tion is controlled by the magnitude of the changes in the 
decision set and the number of iterations required to 
re-allocate the resources throughout all the stages. The 
user must make a value judgement on the importance of the 
optimal solution versus the additional computer costs. 
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In conclusion, the emphasis of this paper is to create 
an algorithm to determine the optimal solution in a quick 
and easy manner. It is evident that a combination of the 
Marginal Analysis and Funnel programs is an efficient 
approach for solving this type of problem. Using the 
Marginal Analysis program to generate a best estimate of the 
solution, the ENTES function with the Funnel program can 
then be applied to confirm the optimal solution or to create 
an input data file that DP6 may use to find the optimal 
solution. The provisioning problems faced by NAVSUP, as 
described in Chap 1, may be solved with this method if they 
can be formatted in the recursive DP construct. The program 
D?6 in particular solves only the minimization of MSRT with 
a constraint on cost. Variations on this program could be 
devised to solve additional objectives like supply material 
availability (SMA ) , time weighted units short(TWUS) , etc., by 
changing the computation procedures in the subroutines. Ihe 
constraint however, is currently restricted to the linear 
form, 

£ Wi*Si<=B i=1,...,N (5.1) 

where Wi is resources required per unit. Si is the number of 
units of item i, and B is the total resource available. 
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4E£INDix a 
DP6 CODE 



This appendix contains the computer listing of the 
program and subroutines for the funneling procedure, 
of the subroutines contain code for objectives 
constraints other than those discussed in this thesis. 



main 

Some 

and 
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